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Abstract. We study dynamics and bifurcations of two-dimensional reversible 
maps having non-transversal heteroclinic cycles containing symmetric saddle 
periodic points. We consider one-parameter families of reversible maps unfolding 
generally the initial heteroclinic tangency and prove that there are infinitely 
sequences (cascades) of bifurcations of birth of asymptotically stable and unstable 
as well as elliptic periodic orbits. 



1. Introduction 

Reversible systems have a very special status inside the realm of dynamical systems. 
Usually, they are positioned "between" dissipative and conservative systems. In the 
context of continuous dynamical systems, a reversibility means that the system is 
invariant under the change of time-direction, t i-> — i, and a transformation in the 
spatial variables. In the discrete context, reversibility of a map / (a diffeomorphism) 
means that / and / _1 possess the same dynamics. Notice that the term "the same" can 
have rather different meaning. If / and f^ 1 are smoothly conjugate, i.e., foh = ho / _1 
and h is a just a diffeomorphism, then / is called weekly reversible. However, much 
more interesting types of reversibility appear when h possesses some structures. For 
example, if h is an involution, i.e., h 2 — Id. In this case, the map / is called 
strongly reversible. Since this last case is the most frequent one in the literature 
(probably, beginning with Birkhoff), nowadays strongly reversible maps are simply 
called reversible maps. 

In contrast to conservative and dissipative systems, the study of homoclinic 
bifurcations in reversible systems is not so popular. Even for two-dimensional 
maps, only few results are known and most of them relate to "conservative and 
reversible" maps which form a certain codimension-oo subclass in the class of reversible 
maps. This situation is probably due to the "common belief" that conservative 
and dissipative phenomena of dynamics only exist separately and, thus, there is no 
necessity to study them "all together" . 

However, they actually can appear together in a dynamical system, giving rise to 
the so-called phenomenon of mixed dynamics, which was recently discovered in [T^] 
(see also [HI [TBI US])- The essence of this phenomenon consists in the fact that 

(i) a dynamical system has simultaneously infinitely many hyperbolic periodic orbits 
of all possible types (stable, completely unstable and saddle) , and 
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(ii) these orbits are not separated as a whole, i.e., the closures of sets of orbits of 
different types have nonempty intersections. 

It was shown in [12] that the property of mixed dynamics can be generic, i.e., 
it holds for residual subsets of open regions of systems. In particular, it was also 
proved that such regions (Newhouse regions, in fact) exist near two-dimensional 
diffeomorphisms with non-transversal heteroclinic cycles containing at least two saddle 
periodic points 0\ and O2 such that | J{0\)\ > 1 and |J((92)| < 1, where J(0;) is 
the Jacobian of the Poincare map (the diffcomorphism iterated as many times as the 
period of Oi) at the point Oj, i = 1,2. 

Let us recall that a heteroclinic cycle (contour) is a set consisting of saddle 
hyperbolic periodic orbits 0\, . . . , O n as well as heteroclinic orbits Tij C W u (Oi) n 
W s (Oj)), where at least the orbits and T n> i for i = 1, . . . , n — 1, are included. 

In general, cycles can include also homoclinic orbits T^i C W u (Oi) n W s (Oi). An 
heteroclinic cycle is called non-transversal (or non-rough) if at least one of the pointed 
out intersections W u (Oi) n W s (Oj) is not transverse. 

If a heteroclinic cycle (or a homoclinic orbit) is transverse, then, as is well-known 
after Shilnikov [36] , the set of orbits entirely lying in a small neighbourhood is a locally 
maximal uniformly hyperbolic set. The situation becomes drastically different in the 
non-transversal case. One can say even that the corresponding system is infinitely 
degenerate, since its bifurcations can produce homoclinic tangencies of arbitrary high 
orders and, as a consequence, arbitrary degenerate periodic orbits [TU1 [TBI 12"0] . 

We remind also that systems with homoclinic tangencies are dense in open regions 
(the so-called Newhouse regions) in the space of smooth dynamical systems [26l [27 | [28 ] . 
Moreover, these regions exist near any system with a homoclinic tangency (or a non- 
transversal heteroclinic cycle). Importantly, these regions are present in parameter 
families unfolding generally the initial homoclinic (or heteroclinic) tangency in certain 
open domains of the parameter space in which there are dense values of the parameters 
corresponding to the existence of homoclinic tangencies. Certainly, such domains are 
called again Newhouse (parameter) regions or Newhouse intervals for one-parameter 
families. In general, it should be clear from the context the kind of Newhouse regions 
considered. 

The existence of Newhouse regions near systems with homoclinic tangencies 
was established in [28] for two-dimensional diffeomorphisms, in [11] [29] [35] for the 
general multidimensional case (including parameter families [28 , 11 ) and in [5] for 
area-preserving maps. The existence of Newhouse regions near systems with non- 
transversal heteroclinic cycles follows immediately from these results, since in such 
case homoclinic tangencies appear under arbitrary small perturbations. Moreover, in 
this case also the so-called Newhouse regions with heteroclinic tangencies can exist. 
It was proved in |12] that if the non-transversal heteroclinic cycle is simple, i.e., it 
contains only one non-transversal heteroclinic orbit and the corresponding tangency of 
invariant manifolds is quadratic, then Newhouse intervals with heteroclinic tangencies 
exist in any general one-parameter unfolding. 

The above-mentioned mixed dynamics takes place as a generic phenomenon. 
This is the case, for instance, when the initial heteroclinic cycle is contracting- 
expanding [12j . that is, when it contains contracting and expanding periodic points 
(i.e. with the absolute value of its Jacobian being greater or less than 1). It is worth 
mentioning that contracting-expanding heteroclinic cycles are rather usual among 
reversible maps. An example of such a cycle is shown in Figure QJa). In this example 
the reversible map has two saddle fixed points 0\ and O2 and two heteroclinic orbits 
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Figure 1. Two examples of planar reversible maps with symmetric non- 
transversal (quadratic tangency) heteroclinic cycles. One can think, for simplicity, 
that the involution R in both cases is linear of form x i— ► x,y i— > — y and, thus, 
FixiJ = {y = 0}. 



r 12 C W u {0 1 ) n W s {0 2 ) and T 21 C W u {0 2 ) n W s {Oi) such that Ji(Oi) = 2 and 
-R(r2i) = T 2 i, R(T 12) = Ti 2 . Besides, the orbit T12 is non-transversal, so that 
the manifolds W u {0\) and W s {02) have a quadratic tangency along Ti 2 . Since 
R{0\) = O2, it turns out that their Jacobians verify J(Oi) = J _1 (02)- If J(Oi) 7^ ±1, 
i = 1,2, then the heteroclinic cycle is contracting-expanding. This condition is robust 
and is perfectly compatible with reversibility. 

Certainly, results of |12) can be applied to reversible maps with such heteroclinic 
cycles and so the phenomenon of mixed dynamics becomes very important and generic. 
However, reversible systems are sharply different from general ones by the fact that 
they can possess robust non-hyperbolic symmetric periodic orbits, more precisely, 
elliptic symmetric periodic points. Thus, one realizes that the phenomenon of mixed 
dynamics in the case of two-dimensional reversible maps should be connected with the 
coexistence of infinitely many attracting, repelling, saddle and elliptic periodic orbits. 
The existence of Newhouse regions (intervals) in which this property is generic was 
already established in [23] for the case of reversible two-dimensional maps close to a 
map having a heteroclinic cycle of the type depicted in Figure Ufa) . 

However, it appears to be true that the phenomenon of mixed dynamics is 
universal for reversible (two-dimensional) maps with complicated dynamics when 
symmetric structures (symmetric periodic, homoclinic and heteroclinic orbits) are 
involved. This universality can be formulated as the following 

Reversible Mixed Dynamics Conjecture Two-dimensional reversible maps with 
mixed dynamics are generic (compose residual subsets) in Newhouse regions in which 
there are dense maps with symmetric homoclinic or/and heteroclinic tangencies. 

We will assume, in what follows, that the involution R is not trivial, i.e. it satisfies 



We will say that an object A is symmetric when R{A) — A. To put more emphasis, 
sometimes the notation self-symmetric may be used. By a symmetric couple of objects 



R 2 = Id, dimFixi? = 1. 



(1.1) 
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Ai,A 2 , we will mean two different objects that are symmetric to each other, i.e., 
R(A 1 )=A 2 . 

Then the symmetric homoclinic (heteroclinic) tangencies from the RMD- 
Conjecture can be divided into two main types: 1) there is a non-transversal symmetric 
heteroclinic orbit to a symmetric couple of saddle points, or 2) there is a symmetric 
couple of non-transversal homo/heteroclinic orbits to symmetric saddle points. 

Notice that the heteroclinic quadratic tangency shown in Fig. HJa) relates to the 
type 1), whereas an example of reversible map having a heteroclinic cycle of type 2) 
is shown in Figure QIb) . The latter map has two symmetric saddle fixed points 0\ 
and 2 (R(0\) = Oi, R(0 2 ) = 2 ) and a symmetric couple of heteroclinic orbits 

r 12 c w u {Oi) n w s (o 2 ) and r 21 c w u {0 2 ) n w^(Oi) (r(t 12 ) = r 21 ). 

As we mentioned above, the case of non-transversal heteroclinic cycles of type 1), 
as in Fig.QJa), was studied in the paper [23] where, in fact, the RMD-Conjecture was 
proved for general one-parameter (reversible) unfoldings, under the generic condition 
J{O l ) = J- 1 {0 2 )^l. 

In the first part of this paper, Sections [2] and |3j we state and prove the RMD- 
Conjecture for one-parameter families which unfold generally heteroclinic tangencies 
of type 2), as in Fig. [ljb). We will call this type as reversible maps with a 
symmetric couple of heteroclinic tangencies. We notice that in this case the condition 
J(Oi) = J{0 2 ) = 1 holds always since Oi £ Fix R, showing that generic conditions are 
different in systems with different types. The generic condition that we are going to 
assume for systems of type 2) is denoted by condition [C] in Theorem[TJ This condition 
(see (|2.ip and comments to it) amounts to say that the (global) map defined near a 
heteroclinic point is neither a uniform contraction (expansion) nor a conservative map. 
More precisely, it has to have a non-constant Jacobian in those local coordinates (near 
0\ and Oi) in which the saddle maps are, a priori, area-preserving. In particular, such 
local coordinates are given by Lemma [2] in which the normal form of the first order 
for a saddle map is derived^ 

Moreover, as we show, symmetry breaking bifurcations have also another nature, 
in comparison with [23J . We find a two-step "fold=> pitch- fork" scenario of bifurcations 
in the first-return maps leading to the appearance of non-conservative fixed points 
which can be either attracting and repelling or saddle with Jacobian greater and less 
than 1 (see Theorem [T] and Figured]). Notice that in the case of heteroclinic cycles 
of type 1), such non-conservative points appear just under fold bifurcations [23] in 
a symmetric couple of first-return maps, whereas elliptic points appear in symmetric 
first-return maps. 

The second part of this work, Section [4] has a more applied character. We show 
in Subsection 14. 1 1 that reversible two-dimensional maps with a priori non-conservative 
orbit behaviour can be obtained as certain periodic perturbations of two-dimensional 
conservative flows of form x = y, y = F(y) . We require that these perturbations 
include explicitly the "friction term" x and preserve only reversible properties of the 
initial flow (for example, they keep the perturbed systems to be invariant under the 
change x v- > x, y t— > —y,t > —t). In this way, in particular, we can obtain the 
reversible maps of type 2), see Figure [2] As a concrete example, we consider, in 
Subsection 14. 1[ the periodically perturbed Duffing equation. In Subsection 14.21 we 

I However, the property of a symmetric saddle periodic point to be a priori area-preserving is 
more delicate. It is well-known, see e.g. |4], that a symmetric reversible saddle map is "almost 
conservative", i.e. its analytical normal form is exactly conservative, and its C°° formal normal form 
(up to "flat terms") is conservative.) 
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consider an example of a reversible map from |30) denned on the two-dimensional 
torus which shows a visible non-conservative orbit behaviour. This map describes 
the dynamics of three coupled simple rotators with small symmetric couplings. The 
couplings are chosen here in such a way that they preserve the reversibility of the 
initial uncoupled three simple rotators. 

In the third part of the paper, at Section [5l we consider a series of problems 
related to the representation of reversible two-dimensional maps in the so-called cross- 
form. It is well-known that such type of cross-forms of maps are very convenient for 
studying hyperbolic properties of systems with homoclinic orbits, both transversal j36j 
and non-transversal [6l O [21] j§| We show that the cross-forms and the corresponding 
cross-coordinates are natural for reversible maps, since they allow expressing many 
reversible structures explicitly and simplify analytical treatment. In particular, in 
the corresponding local cross-coordinates near a symmetric saddle periodic point, the 
normal form of the saddle map becomes very simple. And last (but not least), Section[5] 
contains the proofs of the lemmas needed in the proof of the main theorem [T] 




Figure 2. Examples of reversible maps with non-transversal heteroclinic cycles in 
the cases of (a) "inner tangency"; (b) "outer tangency". These diffeomorphisms 
can be constructed as the Poincare maps for periodically perturbed conservative 
planar systems (left) when the corresponding perturbations preserve the 
reversibility (for example, with respect to the change x — ► x, y — > — y, t — > — t, as 
in the system x = y, y = x — x 3 + ex cost). Notice that the resulting reversible 
map cannot be area-preserving, in general 



§ Since L.P. Shilnikov was the first author introducing such forms and coordinates into dynamical 
systems, they are often referred to as "Shilnikov cross- form" and "Shilnikov coordinates". 
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1.1. Out of the general rule: a collection of reversible maps with codimension one 
homoclinic and heteroclinic tangencies 

It is important to notice that there are many other cases of reversible maps with 
homoclinic and heteroclinic tangencies for which one needs to prove the RMD- 
Conjecture in the framework of one-parameter general families. In Figure [3] we collect 
some simple examples of such maps. They differ by the type of fixed points and 
tangencies: homoclinic or heteroclinic, quadratic or cubicd etc. However, it turns out 
to be more important that all this kind of maps could be separated into two groups: 
a first one including those maps which are a priori non-conservative and a second one 
with those maps where this non-conservativity is, in some sense, hidden. 




Figure 3. Some examples of reversible maps with homoclinic and heteroclinic 
tangencies 

A map as in Figure QJa) (the first one in Figure [3]) belongs to the first group, 
since the condition 7(01,2) 7^ 1 destroys certainly the conservative character. Indeed, 

|| The existence of symmetric cubic homoclinic or heteroclinic tangencies is a codimension one 
bifurcation phenomenon in the class of reversible maps. Therefore, these cubic tangencies should 
be also considered as the main ones jointly with the pointed out quadratic ones. 
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under splitting such heteroclinic cycle, homoclinic tangencies appear both to saddles 
with Jacobian greater and less than 1 and, thus, attracting and repelling periodic 
points can be born [5]. By this principle, all the other cases where a symmetric couple 
of fixed (periodic) points are involved, can be referred to as the class of a priori non- 
conservative maps (e.g., all maps of the first, third and fourth rows of Figure [3]). For 
maps of this type, the problem of finding symmetric periodic orbits (i.e., elliptic ones) 
has to be considered as very important. 

The maps at the second row of Figure [3] have only symmetric fixed points. 
Evidently, the first, third and fourth maps can be assigned to maps with "hidden 
non-conservativity" , since it is not clear, in advance, the existence of bifurcation 
mechanisms leading to the appearance of attracting and repelling periodic orbits. 
It is not the case of the third map (at the second row) which has a symmetric couple 
of (quadratic) homoclinic tangencies to the same symmetric fixed saddle point. One 
can assume, without loss of reversibility, that the map near a homoclinic point is not 
conservative (the Jacobian is greater or less than 1). Then, clearly, stable (unstable) 
periodic orbits can be born under such homoclinic bifurcations. Thus, the main 
problem here is to prove the appearance of elliptic periodic orbits and, as a first 
step, to do it in the one-parameter setting. 

1.2. A short description of the main results 

We describe briefly the central ideas underlying our main results (Theorems [T] and [2]). 
Let assume fo to be a _R-reversible map of type as in Figure [TJb), that is, having 
two symmetric saddle fixed points 0\ and O2 (i.e. R(Oi) = Oi, i = 1,2) and 
two asymmetric non-transversal heteroclinic orbits Ti2 C W u (0±) n W S (C>2) and 
r 2 i C W u {0 2 ) n W s {Ox), satisfying that i?(r 12 ) = T 21 . Let us consider any 
general one-parameter unfolding of fo with /1 being the parameter splitting the initial 
heteroclinic tangency. The main goal of this paper is to show that under general 
hypotheses, any of these unfoldings undergoes infinitely many (in fact, a cascade) of 
symmetry-breaking bifurcations of single-round periodic orbits. Such bifurcations in 
these first-return maps (defined near some point of the heteroclinic tangency) follow 
from the following scenario: 

fold bifurcation =>■ pitch-fork bifurcation 

Under fold-bifurcation, a symmetric parabolic fixed point appears which falls 
afterwards into two symmetric saddle and elliptic fixed points. Concerning (reversible) 
pitch-fork bifurcations, they can be of two classes depending essentially on the type 
of initial heteroclinic cycle they exhibit (see Figure [2] and BJ : in case (a) we say that 
fo has a heteroclinic cycle of "inner tangency" while in case (b) we say that it is of 
"outer tangency" . In the first case, "inner tangency" , from the symmetric elliptic fixed 
point three fixed points are born: a symmetric saddle and asymmetric sink and source. 
In the second case, "outer tangency" , under a pitch-fork bifurcation, the symmetric 
saddle point becomes a symmetric elliptic point and two fixed asymmetric saddle 
points with Jacobian greater and less than 1. Both scenarios are showed in Figure |U 
In the bifurcation diagram of Figure [H related to a conservative approximation of the 
rescaled first-return map, value of c < corresponds to the "inner tangency" case and 

If Notice that in the Lamb-Sten'kin case |23l . see Fig. [TJa), fold bifurcations are directly symmetry 
breaking and lead to the appearance of two pairs of asymmetric periodic orbits: (saddle, sink) and 
(saddle, source). 
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c > corresponds to the "outer tangency" one (situation c = is singular and not 
realisable for our maps; therefore there are no transitions from c < to c > 0). 

We finish this Introduction by describing the structure of the paper. Section [2] 
is devoted to the two above-mentioned main results of the paper, Theorems [1] and O 
Their proof is presented in Section [3] and relies to five lemmas whose proof is deferred 
to Sections [6] and 16.31 Section |4] contains two concrete examples of applications of 
the main results of this paper. 



2. Symmetry breaking bifurcations in the case of reversible maps with 
non-transversal heteroclinic cycles 

Let /o be a C""-smooth, r > 4, two-dimensional map, reversible with respect to an 
involution R satisfying dim Fix(i?) = 1. Let us assume that fo satisfies the following 
two conditions: 

[A] fo has two saddle fixed points 0\ and 2 belonging to the line Fix(-R) and that 
any point Oi has multipliers A,, A" 1 with 0<Aj<l, i = l, 2. 

[B] The invariant manifolds W u (Oi) and W s (0 2 ) have quadratic tangencies at the 
points of some heteroclinic orbit Ti2 and, therefore, by reversibility, the manifolds 
W w (0 2 ) and W s (0\) have quadratic tangencies at the points of a heteroclinic 
orbit T 21 =R(T 12 ). 

Hypotheses [A]-[B] define reversible maps with non-transversal symmetric heteroclinic 
cycles like in FigureQJb). We ask them to satisfy one more condition. Namely, consider 
two points Mi 6 W z " c (Oi) and M 2 € Wf oc (0 2 ) belonging to the same heteroclinic 
orbit r 12 and suppose /<?(Mi) = M 2 for a suitable integer q. Let some smooth local 
coordinates (xi, ?/,•) be chosen near the points Oi in such a way that the local invariant 
manifolds are straightened, i.e., W^ oc (Oi) and W z * c (Oi) have, respectively, equations 
Xi — and y» = 0. Let Ti 2 denote the restriction of the map /q onto a small 
neighbourhood of the point Mi. Then, we assume that 

[C] the Jacobian of Ti 2 is not constant and, moreover, 

Q= 9J ^\m^ (2-1) 

Condition J{T\ 2 ) ^ const is well defined only when certain restrictions on the 
local coordinates hold. One possibility is when these coordinates (xi,y.i) around Oi 
are chosen in such a way that W^JOi) and Wf oc (02) are straightened. However, the 
sign of J{Ti 2 ) depends also on the orientation chosen for the coordinate axes. To be 
precise, we choose these orientations in such a a way that: (i) the y and a;— coordinates 
of the heteroclinic points Mi e Wfc (Oi) and M 2 £ Wf (0 2 ) are positive; (ii) for the 
symmetric points, M{ = i?(Afi) € Wf oc (Oi) and M 2 = R{M 2 ) G ^ c (0 2 ), the x and 
y-coordinates are positive as well. 

Two classes of reversible maps satisfying conditions [A]-[C] can be distinguished: 
those maps with "inner" (heteroclinic) tangency and those with "outer" tangency, 
corresponding to J(Ti 2 ) > and J(Ti 2 ) < 0, respectively. Two examples of such 
diffeomorphisms are shown in Figure [2] Notice that in both cases the global map Ti 2 
is orientable. In the case (a) the axes Xi,yi and x 2 ,y 2 have the same orientation, 
whereas the orientations are different in the case (b). 

Once stated the general conditions for fo, let us embed it into a one-parameter 
family {/^j of reversible maps that unfolds generally at fj, = the initial heteroclinic 
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tangencies at the points of T 12 . Then, without loss of generality, we can take /u, 
as the corresponding splitting parameter. By reversibility, the invariant manifolds 
W u {O x ) and W s {0 2 ) split as W u (0 2 ) and W s {O x ) do when n varies. Therefore, 
since these heteroclinic tangencies are quadratic, only one governing parameter is 
needed to control this splitting. 

Let U be an small enough neighbourhood of the contour C = {Oi, O2, ^12, ^21}- 
It can be represented as the union of two small neighbourhoods (disks) XJ\ and U2 
of the saddles 0\ and 2 and a finite number of small disks containing those points 
of Ti2 and T 2 i which do not belong to U\ and U2 (see Figure [5]). We will focus our 
attention on the bifurcations of the so-called single-round periodic orbits, that is, orbits 
lying entirely in U and having exactly one intersection point with every disk from the 
set U\(Ui U J7 2 )- Any point of a single-round periodic orbit is a fixed point of the 
corresponding first-return map T^ m , that is constructed by orbits of / M with k and m 
iterations (of f^) in U\ and U2, respectively. We will call them single-round periodic 
orbit of type (fc,m). The values of k and m will be always prescribed a priori. The 
first main result is as follows: 

Theorem 1 Let {f^}^ be a one-parameter family of reversible diffeomorphisms that 
unfolds, generally, at fi = the initial heteroclinic tangencies. Assume that fg verifies 
the conditions [A]-[C]. 

Then, at any segment [— e, e] with e > , there are infinitely many intervals 8km 
with border points /4oi'd^ an d Mpf'™' such that 6km — > as k, m — > 00 and the following 



(i) The value fi = /4oici corresponds to a non-degenerate conservative fold 
bifurcation and, thus, the diffeomorphism f ^ has at // € 5km two symmetric, 
saddle and elliptic, single-round periodic orbits of type (k,m). 

(ii) The value fi = fJ^ f ' m ^ corresponds to a symmetric (and non-degenerate if condition 
[C] holds ) pitch- fork bifurcation depending on the type of fo : 

(ii)a In the case of "inner" tangency, single-round asymmetric attracting and 
repelling periodic orbits of type (k,m) are born and, moreover, these orbits 
undergo simultaneously non-degenerate period doubling bifurcations at the 
value (i = Mpd" 1 ' ( w here Mpd" 1 ' as k,m ^ co). 

(ii)b For the "outer" tangency, two single-round saddle periodic orbits of type 
(k, m) with Jacobian greater and less than 1, respectively, are born. Moreover, 
they do not bifurcate any more (at least for < e). 

We refer the reader to Figure Q] for an illustration of this theorem. 

Theorem [T] and its counterpart result in [23] show that the appearance of non- 
conservative periodic orbits under global bifurcations can be consider as a certain 
generic property of two-dimensional reversible maps. 

Briefly, the method we use - based on a rescaling technique - will allow us to prove 
that the first-return map Tfe m can be written asymptotically close (as k, m — > 00) to 
an area-preserving map of the form: 



except 2 = 0. The region c < will stand for the "inner" tangency case and c > for 



holds: 




(2.2) 



in which the coordinates (x,y) and the parameters (M,c) can take arbitrary values 
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Figure 4. Wc mark, by shading, some intervals corresponding to the existence 
of two asymmetric single-round periodic orbits. 



the "outer" one. Its bifurcation diagram is showed in Figure [8] The map (|2.2p is, in 
fact, the product of two Henon maps with Jacobian — c and — c _1 , (sec equations p.4[l ). 
Thus, we can state (see also [ID]) that map (|2.2p has a complicated dynamics in the 
corresponding parameter intervals. The latter means, in particular, that all fixed 
points become saddles and all of them have homoclinic and heteroclinic intersections 
for all values of the parameter (i including (quadratic) tangencies for dense subsets - 
Newhouse phenomenon. 

An analogous "homoclinic tangle" can be observed for the first-return map T/. m 
(see Lemma [3]). However, although the map (|2.2[) is reversible and conservative (its 
Jacobian is identically 1), the original first- return map T^ m is also reversible but not 
conservative in general (see Lemma|4]). Precisely, it will be shown that in some regions 
of the space of parameters (c, M) the map Tfe TO possesses chaotic dynamics and has 
four saddle fixed points, two of them symmetric conservative and a symmetric couple 
of fixed points (that is, symmetric one to each other and with Jacobian greater and 
less than 1, respectively). According to [T2] [23] [20] , the following result holds: 

Theorem 2 Let{f^} be the one-parameter family of reversible maps from Theorem^ 
Then, in any segment [—£, e] of values of fi, there are Newhouse intervals with 
mixed dynamics connected with an abundance of attracting, repelling and elliptic 
periodic orbits. This is, values of parameters corresponding to maps f ^ exhibiting 
simultaneously infinitely many periodic orbits of all these types form a residual set ( of 
second category) in these intervals. 
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3. Proof of Theorem [T] 



3.1. Preliminary geometric and analytic constructions 

To ease the reading all the proofs of the lemmas of this section have been deferred to 
Sections El [6] and 16.31 

Let us consider first the map fo and let M-f G Ux, M% e C/ 2 be a pair of points 
of the orbit T 12 and G U 2 , G U\ be a pair of points of r 21 . Consider 

11+ C Ui and IT" C C/j small neighbourhoods of the heteroclinic points Mf and M~ 
(see Figure [5]). Let us assume that (i) the heteroclinic points are symmetric under 
the involution R , i.e. Mf = R(Mf) and Mf = R(M£), and (m) they are the "last" 
points on U\ and U 2 , that is, fo{Mf) £ U t (and, thus, / _1 (Af+) £ U t ). Let <? be 
such an positive integer that = Jq(M^) (and, thus, — /o(M^~)). 




Figure 5. Schematic actions of the local (Tg, and TqIJ) and global (T12 and T21) 
maps in the neighbourhood U of the contour C = {0\ , O2, I\2, T2i}. 

Consider now the map / M . Denote T oi = f tl \ (J , i — 1,2. The maps T i and T 02 
are called t/ie local maps. We introduce also the so-called global maps T 12 and T 21 by 
the following relations: T 12 = /« : II " -> 11+ and T 12 = /« : il" -> IT+ (see Figure [5J. 
Then the first-return map T]~ m : II + 1— > II + is defined by the following composition of 
maps and neighbourhoods: 

n+ ^ nj n+ ^ n 2 - +H n+ (3.1) 

+ One can always take all neighbourhoods Ui,U~ ,Tlf ,i = 1,2 to be also R-symmetric (that is 
R{U i ) = U i ,R(uf)=nf). 
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Denote local coordinates on Hf and II i as (xoi,yoi) and (xu,yu), respectively. Then 
the chain (|3.ip can be represented (in coordinates) as 

(a;oi,yoi) {xiuVii) ^ (^02,2/02) (2:12,2/12) ^ (£01,2/01)- 

As usually, we need such local coordinates on U\ and Z7 2 in which the maps Tn 
and T02 have their simplest form. We can not assume the maps To; are linear, since 
by condition [A], only (^-linearisation is possible here. Therefore, we consider such 
C r ~ ^coordinates in which the local maps have the so-called main normal form or 
normal form of the first order. 

Lemma 1 (Main normal form of a saddle map) Let a C -smooth map T be 
reversible with dimTz:r(To) = 1. Suppose that T has a saddle fixed (periodic) point O 
belonging to the line Fix{Tq) and having multipliers X and X" 1 , with |A| < 1. Then 
there exist C r ~ 1 -smooth local coordinates near O in which the map Tq (or Tq, where 
n is the period of O) can be written in the following form: 

T _ x = Xx(l + h 1 (x,y)xy) _ _ 

°' V = X~ 1 y(l + h 2 (x,y)xy), 

where h\(0) = — /i2(0). The map i3. 2]) is reversible with respect to the standard linear 
involution (x,y) n- (y,x). 

When proving this lemma one can deduce more "descriptive" properties of this 
map. Namely, one can show that it could be written in the so-called cross-form (see 
Section l5TTj) as follows: 

Tq . x = Xx + h(x, y)x 2 y, _ ^ 

y = Xy + h(y,x)xy 2 . 

Notice that when To is linear, i.e. Tq : x — Xx,y = X~ 1 y, one can easyly 
find formulas for its iterates Tj , j € Z. Namely, one can write it either as 
XjX J Xo,yj = X~iyo or, in cross-form, as Xj — X 3 Xo,yo = X^yj. If To is non-linear, 
then its cross-form equations exist too. In particular, the following result holds: 

Lemma 2 (Iterations of the local map) Let To be a saddle map written in the 
main normal form \3.'J\) (or \3. 3\) ) in a small neighbourhood V of O. Let us consider 
points (x ,yo), • • •) (xj,yj) from V such that (xi +ll yi +1 ) = T (xi,yi), I = 0, - 1. 
Then one has 

Xj = X J x a (l+jX j hj(x ,yjj) , (3.4) 

Vo = ^yj ('i- + jX 3 hj(y j ,x )) , 

where the functions hj(yj, Xq) are uniformly bounded with respect to j as well all their 
derivatives up to order r — 2. 

Remark 1 (a) Both lemmas Q] and [H remain true if To depends on parameters. 
Moreover, if the initial Tq is C r with respect to coordinates and parameters, then 
the normal form is C r ~ 1 with respect to coordinates and C r ~ 2 with respect 

to parameters (see \20l). 

(b) It follows from Lemma[]]that the involution L(x,y) = {y,x) is very convenient for 
the construction of symmetric saddle maps. Moreover, this involution is (locally) 
smoothly equivalent to any other involution R with dim Fix(R) = 1 (see 124-1). 
Thus, our assumption on a concrete form of the involution for the maps f ^ does 
not lead to loss of generality. 
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(c) Similar results related to finite-smooth normal forms of saddle maps were 
established in JS [\ 14\\l^\\20f for general, near-conservative and conservative maps. 
In this paper we, in fact, modify the corresponding proofs adapting them to the 
reversible case. 



3.2. Construction of the local and global maps 

By Lemma [5J we can choose in U\ and U 2 local coordinates (xi,yi) and (£2,1/2), 
respectively, such that the maps Tqi and T02 take the following form: 



101 



xi = Xxxx + h\(xi,yi)x\yi, y x = A x 1 y 1 + hl(xi,yi)xiyf, 



and 



Try 



x 2 = X2X2 + h\(x 2 ,y 2 )x\y 2 , 2/2 = X 2 1 y 2 + h\(x 2 ,y 2 )x 2 y\. 

Furthermore, in these coordinates, the local stable and unstable invariant manifolds 
of both points 0\ and 2 are straightened: Xi = is the equation of Wj^ c (Oi) and 
j/i = is the equation of Wf^Oi), i = 1, 2. Then, we can write the (x, y)-coordinates 
of the chosen heteroclinic points as follows: M± (x^ ,0), M^(0, yj"), M 2 + (xJ,0) and 
M 2 (0, y 2 ). Besides, because of the reversibility, we have that 

xf=y^=at, x^=y 2 =a* 2 (3.5) 

We assume that T 0i (IL?) fMlf = and T^ x (n,r ) nlTr = 0, i = 1, 2. Then the domain 
of definition of the successor map from 11+ into IT~ under iterations of Toi consists 
of infinitely many non-intersecting strips cr? 1 which belong to Hf and accumulate at 
Wi oc (Oi) n n+ as j — > 00. Analogously, the range of this map consists of infinitely 
many strips aj l = T^a® 1 ) belonging to IT~ and accumulating at W^ c (Oi) n FTr as 
j —> 00 (see Figure [6]). 



17(11,) 




Figure 6. A geometry of creation of both domains of definition <r? C 11+ (a) 
and domains of the range a\ C fl — (b) for the maps Ti : fl+ — > IT - . 



It follows from Lemma [5] that the map 
following form (for large enough values of k) 

xn = Afxox(l + kX\h\{xm,yii)) 
Voi = Ajyn(l + kX\h\{yu,XQi)) 



en} i-> crl 1 can be written in the 



rpk 

J 01 



(3.6) 
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and an analogous formula takes place for the map T^ 2 l '■ "m i-> er 



.12. 



T, 



02 



x 12 = \™x 02 (l + mX 2 n h 2 n (x 02 ,yi 2 )) 



2 

2/02 = A^y 12 (l + m\™h 2 m {y 12 ,x 02 )) 
We write now the global map T\ 2 in the following form 



12 



£02 -x% = F X2 {x\\,y\\ -yf = axn + %n 
£02(1/11 - jyf) 2 + ¥>i(a;ii,2m,M)> 

?yo2 = Gi 2 (a:ii,2/ii - yr>A*) = M + c^H + d(yn - 



/ii«ii(yn - ?/i ) + /o3(yn - Vi 

(P2(xn,yn,n), 

where Fi 2 (0) = G 12 (0) = since T 12 (A/~) = Af+ at /x = and 
fi = 0(\y n -Vi\ 3 ) +x n O(\\(xn,yn -Vi)\\), 



-Vi )+ 
Vi?+ (3-7) 



^2 = 0(\ Xll \ 2 ) + 0(\ yil - 2/rl 4 ) + 0(x n (y. 



Vi?) 



Since the curves T i2 (W^ c (Oi) : {^n = 0}) and Wf oc (0 2 ) : {y Q2 — 0} have a quadratic 
tangency at fx = 0, it implies that 



3G 12 (0) 



= 0, 



9 2 G 12 (0) 



= 2d 7^ 0. 



yD 2 



The Jacobian J(Ti 2 ) has, obviously, the following form: 

^ {T\ 2 ) = -be + af iixn + Q (jyii - ) + O (x 2 n + (y n 
where 

Q = 2ad-bf 11 -2cl 02 . (3.8) 
Now condition [C] can be formulated more precisely. Namely, we require that0 
dJ(T 12 ) 



dyn 



¥= 



(3:11=0,^11=^! ,M=0) 

Concerning the global map T 21 , we cannot write it now in an arbitrary form. 
The point is that after written a formula for the map T 12 it is necessary to use the 
reversibility relations to get the one associated to it: 

21— Kl 12 ti , ±i 2 — It ±21 tX 

for constructing T 21 . Then, by (|3.7I) . we obtain that the map T 2 ~^ : n^{(x i, J/01)} ^ 
H2 {(xi 2 ,yi 2 )} must be written as follows 



T 21 < 



£12 = Gi 2 (y i,x i - Vi,^) = 

H + q/ i + d(x i - t/f ) 2 + /ii2/oi(^oi 
/i3(a;oi - yf) 2 + ip 2 (yoi,x i,n), 

yu - x\ = F i2 (yox,xoi-yi, a*) 



Vi 



(3.9) 



Relation 
implicit form. 



ay i +b(xi ) i - y 1 ) + ^02(^01 - Vi ) 2 + <Pl{V0U ^01, A*) 
allows to define the map T 21 : n^Ka;^, j/i 2 )} *— > n^{(a;oi, 2/oi)} ; but in 



* Notice that the analogous phenomenon (of influence of high order terms on dynamics) was 
discovered in 1 151 1171 when studying bifurcations homoclinic tangencies to a saddle fixed point with 
the unit Jacobian. 
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3.3. Construction of the first-return maps T km and the Reseating Lemma 

Now, using relations (|3.6[l - (|3.9[l . we can construct the first- return map T km = 
T 2 iT^T 12 Tq 1 defined on the strip er^ 1 c 11+ . Recall that any fixed point of T km 
corresponds to a single-round periodic orbit of type (k, to) of period (k + m + 2q). 
However, we do not state the problem of studying the maps T km for all large k and 
m. We suppose k and m are large enough integers such that 

A} ~ A£\ (3.10) 

In other words, both values of X k X 2 ,n and X^ k X 2 n are uniformly separated from and 
oo as k, m — > oo. Then the following result holds. 

Lemma 3 (The rescaling lemma) Let the map fo satisfy conditions [A]-[B] and 
fn be a general unfolding in the class of reversible maps. Suppose k and m are large 
enough integer numbers satisfying relation 13.10\) . Then one can introduce coordinates 
(called "reseated coordinates") in such a way that the first-return map Tkm takes the 
form 

M + cy + dx 2 + fuXfxy + f 03 X k x 3 = 

= bX' 2 n X^ k y + aX^x + laiXfy 2 + 0{kX\ k ), 

(3.11) 

M + cx + dy 2 + f\\X\xy + f 03 X k y 3 = 

bX 2 n X^ k x + aX^y + l 2X 2 n x 2 + 0{kXf), 

where 

M = Xf k (/x + cAfaKl + ...)- X r 2 n a* 2 (l + ...)) (3.12) 

and ". . ." stands for some coefficients tending to zero as k,m — > oo. Notice that the 
domain of definition of the new coordinates x ~ AJ~ (xoi — a*), y ~ A^" fc (j/n — a*) 
and parameter M cover all finite values as k,m —¥ oo. 

3.4- On bifurcations of fixed points of the first-return maps T km 

We study bifurcations in the first-return map T km using its rescaled form (|3 . 1 1 [) . If 
we neglect in (|3.1ip all asymptotically small terms (as k,m —> oo), we obtain the 
following truncated form for T km 

M + cy + dx 2 = (3 km y, (i km x = M + cx + dy 2 , (3.13) 
where /3fe m = 6A^ fc A™. Rescale the coordinates 

fikm i^km 
X j X new 7 y — j ynew- 

d d 
Then map (|3 . 1 3[) is rewritten in the following form, that we denote with H: 

H : M + cy — x 2 = y, x = M + cx- y 2 , (3.14) 

where 

1 

PL. "" " P~m~b- 



M = -&- M , 2= — = I A^A 2 - m (3.15) 



Notice that H depends on two parameters M and c which can take arbitrary values, 
except for c = (according to conditions be ^ and < |Aj| < 0). Thus, two mainly 
different scenarios take place: with c < and c > 0. 
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Observe that the map H can be expressed in the explicit form (I2.2[) . Moreover, it 
can be represented as the superposition H = I-L2 of two quadratic (Henon) maps 

"Hi : { x x = y , yi = M + cx - 

n, r Ml 1 2 

«2 : { #2 = 2/1) 2/2 = — - + -xi + -2/1 • 

c c c 

The Jacobians of these maps are constant and inverse: J{ r H\) — — c and J{J-L%) = 
— c . Therefore, the resulting map H — %2°Hi is a quadratic map with the Jacobian 
equal to 1 and so area-preserving|fl 

Form (I3.14p of H allows to give a rather simple geometric interpretation of the 
bifurcations of fixed points. The coordinates (x, y) of these fixed points must satisfy 
the equations 

2/(1 -c) = M-x 2 , x(l — c)=M — y 2 (3.16) 

Let us hold fixed c and suppose that c ^ 1. Then equations of (|3.16[> define on the 
(x, 2/)-plane two parabolas which are symmetric with respect to the bisectrix y = x. 
Intersection points of the parabolas are also fixed points of H. When M varies the 
parabolas "move" and, as a result, the number of intersection points can change (i.e. 
bifurcations in H occur). See Figure [7] in which the case c < 1 is illustrated: (a) 
the parabolas do not intersect if M < = — |(c — l) 2 ; (b) the parabolas are 
touched (quadratically) to the bisectrix and one to other, if M = M*; (c) they have 
two (symmetric) intersection points if M* < M < M% = |(c — l) 2 ; (d) they have 
a cubic (symmetric) tangency when M = and, finally, (e) the parabolas have 
four intersection points (two symmetric points and a symmetric couple of points) if 
M > M|. An analogous picture takes place for the case c > 1 (the parabolas have 
their branches in the opposite directions). The case c = 1 is very special. Here, the 
equation (|3.16[) takes the form = M — x 2 , = M — y 2 and, thus, a certain "0 — 4"- 
bifurcation occurs at M = 0: the map H has no fixed points for M < and 4 fixed 
points appear immediately when M becomes positive. 

More details concerning bifurcations of fixed points of H are illustrated in 
Figure [8] where principal elements of the bifurcation diagram on the (c, M)-planc 
are represented. Notice that the line c = is singular and, therefore, there are no 
transitions between the half-planes c < and c > 0. In particular, this means that 
bifurcation curves must "terminate" on the line c = (two such terminated points 
are denoted in Figure [8] as black stars). Besides, three types of bifurcation curves are 
represented in the figure: fold (F), period-doubling (PD) and pitch- fork (PF). 

The curves Fi and F2 having the same equation 

M = -~(c -l) 2 (3.17) 

but with c < and c > 0, respectively, relate to a conservative fold-bifurcation. If 
M < -i(c-l) 2 , i.e. (c,M) £ Ii Ul r , the map H has no fixed points; i£M > -\(c-l) 2 , 
the map H has two symmetric fixed points P + = (p+,p+) and P~ = (p_,p_), where 



c-l± J c-l) 2 +4M . 1 
P ± M = — 4 (c- I) 2 (3.18) 

(I Notice that a non-trivial technique proposed in 1401 based on considering superpositions of Henon- 
like maps, allows to deduce a series of quite delicious generic properties demonstrating richness of 
chaos in non-hyperbolic area-preserving maps. 
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(a) M < 



(b) M = M 1 * 



(c) M * <M<M 





(d) M=M * 



(e) M>M 



Figure 7. A geometric interpretation for a structure of fixed points of map (13 - 141) 
as intersection points of two symmetric parabolas. 



Notice that the point Q* = (c = 1,M = 0) 6 F2 corresponds to a degenerate fold- 
bifurcation: simultaneously four hxed points, two symmetric and a symmetric couple, 
are born at the transition I r — > V r . 

In the case c < 0, the point p_ is always a saddle and it does not bifurcate any 
more. This is not the case of the point p+ which can undergo both period-doubling 
and pitch-fork bifurcations. The period-doubling bifurcation curves 



PD\p, 
PD 2 (p, 



M 



M = 



1 



l--(c-iy, c<0; 



are represented in Figure [5] by "grey arrows" which indicate directions of birth of 
period-2 points. The curve 

PF 1 ^): M=|(5-l) 2 , c<0 

relates to the pitch-fork bifurcation: when crossing this curve (in the direction 
V; — > VI;) the point p+ becomes saddle and two asymmetric elliptic fixed points 
P3 and pi are born in its neighbourhood. The point p + does not bifurcate any more, 
whereas, the points ps and pi undergo simultaneously period-doubling bifurcation at 
crossing the curve 

(1 - 3c) (3- 5) 



PD(p 3A ): M 



c < 



Further variation of parameters in the domain VII;, will lead to a cascade of 
(conservative) period-doubling bifurcations of asymmetric periodic points. 
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Figure 8. Elements of the bifurcation diagram for the map H. 



As it can be seen in Figure the character of the bifurcations in the case c > 
is different to the one for c < 0. In this case, c > 0, both symmetric fixed points 
p+ and p- undergo pitch-fork and period-doubling bifurcations. The corresponding 
bifurcation curves are 

PF 2 {p+): M=^(c-1) 2 , 0<c<l; 

PD 3 (p+): M = l-l( c -lf, 5>0; 

for the point p+ and 

PF(p_): A} = 5(c-1) 2 , E> 1j 

PD(p _ ): M = <° +1 >f-", £ >0 <3 ' M) 

for the point p_. Notice that, when crossing the curve PF 2 (p + )UPF(p_), a symmetric 
couple of saddle fixed points p$ and P4 are born and they do not bifurcate any 
more (for values of parameters c > and M in the domain lying above the curve 
PF 2 (p + ) U PF(pJ)). One can expect, however, that bifurcations of the symmetric 
fixed points p+ and p- give rise to cascades of period-doubling bifurcations. 
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It should be noted that, despite the reversibility, the pitch- fork and period- 
doubling bifurcations can have in T^ m a different character in comparison with the 
one of the truncated map H. Of course, the pitch-fork bifurcation in T^ m leads again 
to the appearance of two non-symmetric fixed points and P4 but these points can 
be non-conservative. In fact, this is a general property (it holds for open and dense 
set of systems) that the following lemma shows. 

Lemma 4 (Non-conservative fixed points) The non- symmetric fixed points P3 
and P4 of map p.llj) . with k and m satisfying p.lOj) . have Jacobian J ns and J„}, 
respectively, with 

J™ = 1 + Q(T? *~ r) A^ + o(A^), (3.20) 
be 

where Q is the coefficient given by formula (|3 .8[) and £* and rj* are, respectively, the 
x- and y-coordinate of the fixed point. 

Due to the reversibility, the fold bifurcation in the first-return map T^ m has the 
same character as in the truncated map H and leads, therefore, to the appearance 
of two symmetric fixed points, p+ and p_, saddle and elliptic ones. Concerning the 
symmetric elliptic fixed points, we have the following result. 

Lemma 5 (Symmetric elliptic fixed points) The point p+ (resp. the point p_) 
is generic elliptic, that is, it is KAM-stable, for open and dense sets of values of the 
parameters (c, M) in the domains Hi U V; and IV r U V r U VI r ( resp. in the domain 
II r UV r U VHIr). 

The proofs of lemmas [3HS] are given in Section 16.31 



3.5. End of the proof of Theorem^ 

If k and m are large enough and having in mind (|3.10[) , (|3.12|) and (|3.15[) , the following 
relation between the parameters p and M holds: 

p = y 2 n a* 2 (l + p\) - cXlal(l + p 2 k ) - ^MA| m (l + pi), 

where p\,i — 1,2,3, are some small coefficients (pi — > as k — > oo). Using formulas 
(|3 . 1 T[) (|3 . 19[) for the bifurcation curves of the truncated map (|3 . 14[) . asymptotically 
close to the (rescaled) first-return map p.llj) . we find the following expressions for 
the bifurcation values of p at the statement of Theorem [1] a value p — p^[^ , which 
corresponds to the fold bifurcation in Tfc m , 

^L m) = A?o5(l + p\) c\\a\{l + pi) + 

^(b-cX^fxrd + pl) 

and a value p = p^J m \ associated to the pitch-fork bifurcation in T^ m (which is not 
conservative if we have in mind Lemma U]) , 

4T } = A?o5(l + p\) - cXlat(l + pi) - 

^{b-cx^fxra+pD. 

These considerations imply Theorem [TJ 
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4. On applied reversible maps with mixed dynamics 

In this section we present two concrete examples of reversible systems where Theorem[T] 
applies and exhibiting, therefore, mixed dynamics: periodically perturbed Duffing 
equation and the Pikovsky-Topaj model |30) for coupled rotators. 



4-1- A periodic perturbation of the Duffing Equation 
Let us consider the following system 

x = y, 



y = — x + x 3 + e (a + (3y s'mujt) , (^-l) 

where a, j3, w € R and e is an small perturbation parameter. The unperturbed system, 
for e = 0, corresponds to the so-called Duffing system (also called Anti-Duffing for 
several authors). It is Hamiltonian, with 

ip- x^ 1 I 

H{x,y)= y -+V(x), V {x) = ----- 

and is (time)-reversible with respect the following linear involutions R(x, y) = (x, —y) 
and S(x,y) = (—x,y). This system has three singular points: one elliptic at (0,0) 
and two saddles at (±1,0). Moreover, these two points are connected through two 
(symmetric) heteroclinic orbits F h . 

The perturbed system, for e ^ 0, is still i?-reversible but, in principle, non 
necessarily Hamiltonian. This is a particular case of a more general family of ir- 
reversible perturbations 

& = y, 

is = -x + x 3 +eg(x,y,t), 

satisfying that g(x,—y,—t) = g(x,y,t). It is well known that two (symmetric) 
hyperbolic periodic orbits 7^ appear close to the saddle points (±1,0). Let us denote 
by W u,s (^) their corresponding unstable and stable invariant manifolds, respectively. 
Generically these invariant manifolds will intersect each other transversally and will 
remain close to the unperturbed heteroclinic connection. The first order in e associated 
to their splitting, will be given by the well-known Poincare-Melnikov-Arnol'd function 

/+00 
L F G(T h (t)) dt, 
-00 

where F(x,y) = (y,—x + x 3 ), G(x,y) — e(0, a + fHysmuit), Th(t) is any of both 
unperturbed heteroclinic connections rr(t) and Lp(G) = (DF)G stands for the Lie 
derivative of G with respect to F. Simple zeroes of M(to) provide tangent intersections 
between the invariant manifolds W u ' s '(7| t ) ■ This systems constitutes a good candidate 
to apply our results. 

For the computation of M(to) we consider here the positive heteroclinic orbit 
Fh = but, by symmetry, everything applies exactly for . Thus, 

T h (t) = (x h (t),y h (t)) = (x h (t),x h (t)) = ftanh-^=, -^=scch 2 --^= j 
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and 

/+oo 
(DF)G\ {Xh(thyh{tht+to) dt= (4.2) 
-oo 

yh{t) (a + —j= ( sech 2 -^= ) sino;(i + to) J dt 



v 2 \ V2 

ea f +oc ,,i , e/3 /" +oc / , 4 i \ . , 

— r= / seen — p at + — / I seen — p I sinw(t + to) ere = 

V2 J-oo v2 2 J_ 00 \ ^2/ 

+ CO) 

Concerning Ji it is straightforward to check that its value is 2. Regarding ^2, it is 
more convenient to compute the integral 

f + OC 



using the method of residues. Indeed, from it we can derive that 

+co/ , 4 f\, , . , 2ttw 2 (lj 2 + 2) 

sech — = sinw(i + to) dt = — . ^ — smwto, 



V2/ 3 sinh^§ 

f + °° ( ,4 * \ , s , 2^w 2 (w 2 + 2) 
loo I 8 71 J C0SW( * + h) = T sinh^ C0S ^° 
and, substituting in (|4.3p . we get 

3sinhp| 

V2 e(aV(uj) + /3sinwt ) , 



w 2 (w 2 + 2) 
provided we dehne 

V2lu 2 (lu 2 + 2) 



3 sinh ^ 

V2 



Therefore, for small values of e we have: (i) if |/3/a| > V(co) then ^"(7") and iy s (7+) 
intersect; (ii) if < V(u>) they do not intersect each other and (Hi) if \/3/a\ = V(lS) 
then M(t ) has zeroes which are double but not triple since dM(t(,)/da = y/2 =/= 0; 
this case leads to quadratic heteroclinic tangencies. 



4-2. On the Pikovsky-Topaj model fSOj/ of coupled rotators 
Let us consider the following system 

ipl = 1 — 2s sin ipi + s sin ip2 

if>2 — 1 — 2e sin ip 2 + s sin ipi + e sin ip 3 (4.4) 
^3 = 1 — 2esinV'3 + e sin -02, 

where ipi £ [0,27r),i = 1,2,3, are cyclic variables. Thus, the phase space of (|4.4[) is 
the 3-dimensional torus T 3 . System (|4.4j) is reversible with respect to the involution 

TZ: Ipl -> 7T - ^3 , V>2 -> 7T - "02 , - 03->7T-V'l- 
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System (|4.4[) was suggested by Pikovsky and Topaj in the paper [3U] as a simple 
model describing the dynamics of 4 coupled elementary rotators. By means of the 
coordinate change 

01-03 01 + ^3 - 7T 01 + 1p 3 - 7T 

e=^— , »?= g , P= g +02-- 

and the change in time dr = dt(2 + ecos(p — rj)) system (|4.4I) is led into 
• 2esin^sin77 
2 + ecos(p — 77) 
1 — ecos(p — rj) — 2ecos£cos77 
? ' ~ 2 + ecos(p-^) ^ ' ^ 

P=l 

Then time-1 Poincare map of system (j4.5[) is also reversible with respect to the same 
involution R : £ — > £, 77 — > —77. 

It was found in [30] that, for small e, system (|4.4p behaves itself as a conservative 
system close to integrable one and several invariant curves could be observed. However, 
when one increases the value of e invariant curves break down and chaos appears 
(which is already noticed, for instance, at e « 0.3). This picture looks to be quite 
similar to the conservative case. However, certain principal differences take place. In 
particular, a "strange behaviour" of the invariant measure is observed. Iterations of the 
initial measure are convergent to some suitable limit. However, the limits t +00 and 
t — > —00 for the same initial measure are different ( numerically observed, for instance, 
for values of e w 0.3). This situation is impossible when the invariant measure is 
absolutely continuous. Therefore, it must be singular and concentrated on "attractors 
and conservators" at t — > +00 or "on repellers and conservators" at t — > —00. Here 
under the term "conservator" we mean the set of self-symmetric non-wandering orbits. 
Moreover, +00- and — oo-invariant measures look like symmetric (with respect to the 
fixed line of the involution) and having non-empty intersection so there are no gaps 
between asymmetric and symmetric parts. This means that "visually" attractors and 
repellers intersect and it is an evidence of mixed dynamics in this model. 

Moreover, a transition from conservative dynamics to non-conservative one can 
be generated by bifurcations of periodic orbits. For small enough e periods of all such 
orbits are large and the corresponding resonance zones are narrow. When increasing 
£, periodic orbits of no too large period appear and dissipative phenomena can become 
observable. For example, the map T under consideration has no points of period 1 
and 2 for e < 0.6 but it has, at e = e* s» 0.445, two period 3 orbits. Notice that 
these orbits are different since map T has the symmetry £ — > 2n — £ that implies the 
appearance of 2 (in fact, an even number) different orbits. Thus, the scenario is the 
following: there is no fixed point for T 3 at e < e*; at e = e* two fixed points with 
double multiplier +1 appear in Fix R one symmetric of each other; at s > e* all 
these orbits fall into four orbits, two symmetric elliptic and two asymmetric saddle. 
Moreover, the latter orbits satisfy that the Jacobian is greater than 1 at one point and 
less than 1 at other point. 

Bifurcations of such type (i.e., having a single point falling into 4 points) are 
not typical in one-parameter families even in the reversible case. Here, general 
bifurcations are met (for symmetric fixed points) of types "0 — > 2" or "1 — > 3" , 
that is, "conservative" fold and "reversible" pitchfork, respectively. The presence of 
a typical bifurcation "0 — > 4" says us about the existence of a certain additional 
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degeneracy in the system. The "clear symmetry" £ — > — £ is not suitable for this role. 
However, system (|4.5[) possesses such a "hidden symmetry" which implies that the 
map Tf p — \_^fp—2ir) i s the second power of some non-orientable map. This peculiarity 
is caused by the fact that the maps T , ( p=w )^( p=27r ) and Ti p= q\^.i p=7 a are conjugate. In 
particular, one can check that 

T(/)=7r)->(/)=27r) = S 1-1 T'(p=0)^(p=7r) S, (4.6) 

through the linear change of coordinates £ — > ir — £, rj — > 77 + 7r, p — > p + ir. Indeed, 
after this coordinate transformation, the right sides of system (|4.5p remain the same, 
but the limits of integration (along orbits of system (14.51) to get the correspondence 
map between sections p = a and p = b) are shifted in tt. Such a property is called 
time-shift symmetry. 

iFrom ([H) it follows that T (p=0) ^( p=27r ) = T (p=0) ^( p=Tr) ST (p=0) ^ (p=x ) 
Since S 2 = Id, one has that S = S^ 1 and, therefore, 

T (p=0 )^(p=27 r ) = (T , ( p=0 )^ (p=T )5) 2 (4.7) 

This means that the map T( p=0 )_>.( p= 27r) considered is the second power of some map. 
Notice that the transformation associated to S is non-orientable and, thus, the map 
T( P= Q)^.( p= T T )S is non-orientable as well and, on its turn, our first-return map T is also 
the second power of some non-orientable map. 

It is straightforward to check that the map T^ p= ^^^ p= ^ is reversible with respect 
to the involution R\{x,y) — (—x,—y) and that the map T( p= o)->( p= „) S is reversible 
under the involution R(x,y) = (x,—y). Thus, the bifurcation of map T 3 at e = e* 
can be treated as a bifurcation of a fixed point with multipliers (+1,-1) in the case 
of a non-orientable map (in fact, the map (T(p = o)->(p=7r) S) 3 ). So, summarising, in our 
case this bifurcation leads to the appearance of two elliptic points of period 2 on Fix R 
and a symmetric couple of saddle fixed points (that is, outside Fix R and symmetric 
one to each other). These saddle fixed points are not conservative. It can be checked 
numerically that the Jacobian of one point is greater than 1 and less than 1 at other 
point. Due to reversibility, the stable and unstable manifolds of saddles pairwisc 
intersect and form a "heteroclinic tangle" zone. This zone is extremely narrow since 
the separatrix splitting is exponentially small. However, moving slightly away from 
the bifurcation moment we can find numerically heteroclinic tangencies and, hence, 
moments of creation of non-transversal heteroclinic cycles. Since the saddles involved 
are not conservative, it follows from |23j the phenomenon of mixed dynamics. 

5. Cross-form type equations for reversible maps. Proof of Lemmas [TJ 
and [2] 

5.1. Cross-form for reversible maps 

As it will be seen along this section, the so-called Shilnikov cross-form variables 
constitute an essential (and natural) tool to deal with reversible maps and a simple 
way to generate them. The first part will be devoted to introduce such variables and 
to present some of its main characteristics. In the second part we apply them to prove 
Lemmas [T] and [21 

We say that a map is in cross-form if it is written as 

( x = h(x,y), 
\ V = h(x,y), 
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On the other hand, let us consider a diffeomorphism F of the plane which is reversible 
with respect to a (in general, non-linear) involution R (R 2 — id, R ^ id), having dim 
Fixi? = 1. Moreover, let assume that the involution R reverses orientation (that is, 
det DR < 0), which is the most common situation in the literature. Our aim is to show 
how reversible maps can be expressed in cross-form type equations and, conversely, 
how reversibility can be derived from this form. 

As an starting point, let us consider the linear set up, that is, when the reversor 
R is the linear involution L : (x, y) n- (y, x). In this case, the following result holds: 

Lemma 6 Any diffeomorphism F : (x, y) t-¥ (x, y) defined, implicitly, by means of 
equations of type 

F ■ < ^ = (5 l) 

1 y = f(y,x) 

is always reversible with respect to L(x,y) = (y,x). 

Proof. Remind that if G is a L-reversible diffeomorphism it must satisfy that 
G o L o G = L or, equivalently, L o G o L = G _1 or (L o G o L) 1 = G. In our 
case we will prove that F, defined by (|5.1|) . verifies the latter relation for G = F and, 
consequently, is L-reversible. To do it, we will use an equivalent expression for the 
inverse of a planar diffeomorphism. Precisely, if 

H f x = hi{x,y), 
\ V = h 2 {x,y), 

the corresponding inverse map H^ 1 : (x, y) i— >• (x, y) can be implicitly written through 
the expression 

H -i . f x = hi(x,y), 
\ y = h 2 (x,y). 

This is clear since (x,y) = H(x,y) implies that (x,y) — H~ 1 (x,y). An algorithmic 
way to get it consists in swapping bars among the variables, that is x ■(-> x and y ■<-> y. 
We apply this procedure to compute formally an expression for (LoFo L)^ 1 and to 
check afterwards that it coincides with F . Let us compute it step by step. First we 
have 

f o l ■. ( * =f }?jy\ 

{ x = f{y,y). 

To apply L onto F o L corresponds to swap x ■(-> y in the precedent expression: 

LoFoL- { y = f( y J^> 
' \ x = f(x,y). 

And finally, to get its inverse we swap bars and no-bars, that is x O x and y ■<-> y. 
Performing this change we obtain 

(LoFoL)- 1 : { V = f f f X ): 
\ x = f(x,y), 

which is exactly the expression for F. So, F given in the form (|5.ip is always ir- 
reversible. 

□ 



This result can be useful to provide suitable local expressions for reversible 
diffeomorphisms in the plane. Thus we have the following lemma. 
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Lemma 7 Let F = (/i,/2) be a planar diffeomorphism, reversible with respect a 
general involution R, C r , r > 1, orientation reversing and with dim Fix R = 1. Let us 
assume the origin (0,0) a fixed point of the involution R, that is (0,0) € FixR. 

Then, if D xx f\ + D yy f2 ^ at (0,0) there exist local coordinates, that we denote 
again by (x,y), in which F admits the following implicit (normal) form 

( x = g(x,y), 

X y = g{v,x). 

This map is reversible with respect to L(x,y) = (y,x). 

In the case of a saddle fixed point, the concrete type of implicit normal form that can 
be obtained is given in equation p.3[) . 
Proof. This will be achieved in two steps: 

(i) First we apply Bochner Theorem [2] which allows us to conjugate, around (0,0), 

our involution R to its linear part DG|( ,o)- 
(ii) Using that the partial derivatives on (0, 0) do not vanish simultaneously, we apply 
Implicit Function Theorem to reach the final form. 

We proceed as follows: 

(i) Notice that if R is a (general) involution and p S Fix R then its linear part DR\ p 
is an involution as well. Indeed, 

Id = R 2 I = D(R 2 )(p) = DR\ R[p) ■ DR\ p = (DR\ p f . 

Bochner Theorem ensures the existence of a C-diffeo ip which conjugates, locally 
around p, R to DR\ P . We include, for completeness, a simple proof of this fact 
given in |33j . From the equality 

DR\ p o(R + DR\ p ) = DR\ p o R + id = 

id + DR\ p oG = (G + DG\ P ) o G 

it follows that DG\ p o(R + DR\ p ) = (R + DR\ p )oR. We define ip = R+DR\ p e C r 
and check that it is a diffeomorphism in a neighbourhood of p: 

Dif}\ p = D(R + DR\ P ) \ p = DR\ p + DR\ p = 2DR\ p 

and so det Dip\ p = 2det DR\ p ^ since R is a diffeomorphism around p. So tjj 
is a C r -diffeomorphism which conjugates R to DR\ p around p. 
Since R is orientation reversing its linear part around p, DR\ p , is also orientation 
reversing. Following |33j for instance, we know that there exists a transformation 
which conjugates DR\ p to the linear involution L(x, y) — (y, x), which will be the 
one we will consider, locally, from now on. 

(ii) Let us assume, for instance, that D yy f ^ at (0,0). Using Implicit Function 
Theorem, we can write from equation y — /2 (#,£/) an expression for y, say 
y = g(y,x), for a suitable function g. Substituting it into the equations defining 
F we get a (locally) equivalent expression for F: 

F . f x = fi(x,y) = fi{x,g(y,x)) =: h(x,y), 

As stated above, we can assume F to be locally conjugated around the origin 
to the linear involution L : (x,y) i— > (y,x). So in that variables (to simplify the 
notation we keep the same name for the variables and the functions involved) 
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it must satisfy that (LoFoL) -1 = F. Applying the procedure introduced in 
Lemma |51 one obtains that 



FoL 



x = h(y,y), 
x = g(y,y)- 



We apply L (that corresponds to swapping x and y, 
LoFoL 



JJ = h(y,x), 
x = 9(x,y) 



and, finally, we swap (x,y) for (x,y), 



{LoFoL 



-l 



V = h(y,x), 



Since it must coincide with F it turns out that h(x,y) — g(x,y) and so 



F 



x = g(x,y), 
y = g(y,x). 



□ 

We present now a counterpart result when the map is given in implicit form. 

Lemma 8 Any map G : (x,y) i— > (x,y) given by 

g(x,y,x,y) = 0, 
g(y,x,y,x) = 0, 

is L-reversible, where L : (x,y) H ^ (y, x). The second equation g(y, x, y, x) = is a 
kind of L- conjugate of the first equation g(x,y, x,y) = 0. 

Proof. It is enough to check that L o G o L = G _1 . To do it we proceed again as in 
Lemma [S] First, remind that an implicit expression for G _1 is always obtained by 
swapping bars for no-bars, that is, (x, y) O (x,y). So 

g{x,y,x,y) = 0, 
g(y,x,y,x) = 0. 

On the other hand wc compute L o G o L. Thus, 

g(y,x,x,y) = 0, 
g(V,x,x,y) = 0, 

and, swapping (x,y) for (y,x), we get 

g(y,x,y, y) = 0, 
g(x,y,x,y) = 0, 

which coincides with G _1 . Therefore the lemma is proved. 

□ 



G- 1 



GoL 



LoGoL 



The following result establishes an interesting relation between polynomial 
reversible and area preserving maps. 

Lemma 9 (|32j) Any Taylor truncation of a planar polynomial diffeomeorphism 
which is reversible with respect to a linear involution is area preserving. In particular, 
this applies to the truncation of a normal form of such diffeomorphisms. 
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Proof. Let z = G(z) a polynomial planar map which is reversible with respect to a 
linear involution S (S 2 = Id, S 7^ Id). This means that S o G o S = G -1 and, in 
particular, that G _1 is also a polynomial. Differentiating the latter expression we get 

SDG\ Sz S = Die- 1 )], = (DGla-i^y 1 
det(^gl,^) = ^-— . 

Using that det^DG^ S) = (det S") 2 det DG\ Sz = dctDG| S2 it follows that 

(detDG| 5l )-(detDG| G -i W ) = l, Vz. (5.2) 

Since G and G" 1 are polynomials and S linear we obtain that det-DGIs^ and 
det_DG|G-i( z ) are polynomials as well. But the product of two polynomials is a 
constant if and only if they are constant, that is, det_DG| 2 = k = constant. Thus, 
from (|5.2|) it follows that k 2 = 1 and, therefore, det_DG| z = ±1, Vz. 

□ 



And last, but not least, we remark another interesting property regarding this 
cross-form type: any polynomial truncation of a reversible diffcomorphism written 
in cross-form type is also in cross-form type and, consequently, it is reversible. This 
means, from Lemma O that this truncation is always area-preserving. 

5. 2. Proof of Lemma [7] 

Let O be a fixed saddle point of a reversible map To. Applying Bochner Theorem [24], 
we can assume the existence of local coordinates around O such that O is located at 
the origin and that the involution R is exactly (x, y) H> (y, x) in these coordinates. 

Let x — v(y) be the equation of the stable manifold. Then, by the i?-reversibility, 
y = v(x) is the equation of the unstable manifold. If \dv/dy\ < 1, we perform the 
transformation x new — x — v(y), y n &w = y — v{x), while, if \dv/dy\ > 1, the change 
is x new — y — v{x), y n ew = x — v(y). After such transformation, which commutes 
with R, the equations of the stable and unstable manifolds become y — and x = 0, 
respectively. Thus, in the corresponding local coordinates, the map can be represented 
in the following form 

x = \x + gi(x,y), y = A _1 y + g 2 (x,y) (5.3) 

where gi(0, y) = 0,g2(x,0) = and g-(0,0) = 7 i = 1,2. It is very convenient to 
rewrite this equation in the so-called cross-form: 

x = Xx + gi(x,y), y = Xy + g 2 (x,y) (5.4) 

Equation (|5.4p comes from (|5.3[) writing y — F(x,y) (which exists due to the Implicit 
Function Theorem) and substituting it into the first equation: x — Ax + 5i(x, F{x, y)). 
The i?-reversibility of (|5.4j) implies that g\{x,y) = g2(y,x) so we can represent map 
(|5.4j) in the form 

x = \x + ipi (x) + ipi (y)x + pi (x, y)x 2 y , g g . 

y = ^y + vi{y) +^i(x)y + pi(y, x)xy 2 1 '' ' 

Perfoming the i?-invariant change of variables 

£ = x + xhi(y), i] = y + yh 1 (x) (5.6) 
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with hi(0) = 0, it turns out the following equation for £: 
£ = x + xhx{y) = 

A£ - xhi (y) + xi/ji (y) + {Xx + ipi (y)x + ipi (x))hi (y) + 

AC + ^i(0 + O(e 2 y) + 

x {-hi{Xy + <fi{y)) + Mv) + + Mv))hi{y)\ ■ 

Since we want the expression in the square brackets to vanish identically, we ask the 
function hi(y) to satisfy the functional equation 

hx(Xy + = M£)(i + A-Vi(y)) + A~Vi(y), (5.7) 

which has solutions hi = hi(u) in the class of C"* -1 -functions. Indeed, we can consider 
(|5.7p as an equation for the strong stable invariant manifold of the following planar 
map 

h = fti(l + A- 1 V'i(u)) + A-Vi(w), 
u = Xu + ip\(u). 

Since < |A| < 1, ip\ (0) = and <fi(0) = y'i(0) = 0, this map has strong stable 
invariant manifold W ss passing through the origin, that is, satisfying an equation 
hi = hi(u) with hi(0) = 0. Therefore, after the i?-invariant change (|5.6I) . the map 
(|5.5p takes the form 

x = Xx + ip 1 (x)+ p 2 {x,y)x 2 y _ g . 

y = H) + vi{y) + P2{y,x)xy 2 

Applying a i?-invariant change of variables of the form 

£ = x + h 2 (x)x, r] = y + h 2 (y)y 

with /i2(0) = 0, the first equation of system (|5.8p can be rewritten, in these new 
coordinates, as follows 

£= X£_+x[-Xh 2 (x) +<pi(x)+ (5.9) 
h 2 (Xx + <p(x))(X + cpi(x))} + 0(£ 2 r/), 

where we have denoted (fi(x) = <f>i(x)x. As we did above for hi, we seek for a function 
h 2 satisfying the following equation 

h 2 {Xx + <p(x)) = (1 + X-^iix))- 1 ^^) - X~ 1 i^i(x)), (5.10) 

which vanishes the expression inside the square brackets in (|5.9[) . As before, equation 
(|5.10|) has solutions h 2 — h 2 (u) in the class of C"" _1 -functions. Again, one can consider 
the expression (|5 . 10[) as an equation for the strong stable invariant manifold associated 
to the following planar map 

h 2 = (1 + X~ 1 {pi(x))~ 1 (h 2 - X-^i(x)), 
u = Xu + (fii(u). 

Having in mind that < |A| < 1 and </?i(0) = ^4(0) = 0, this map admits strong 
stable invariant manifold W ss passing through the origin, i.e., having an equation 
h 2 = h 2 (u) with h 2 (0) = 0. This completes the proof of the Lemma. 
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5.3. Proof of Lemma\^ 

We write the map T in the following form 

x = Xx + h(x,y), y = jy + g(x,y) 
where we assume that 7 = A -1 and 

h(x, y) = x 2 y(p 1 + 0(\x\ + \y\)), g(x, y) = xy 2 {fo + 0(\x\ + \y\)). 
Consider the following operator $ : [(xj, yj)] k =0 >-> [(%!%)] j=o : 



J- 1 ._ _ . 
x 3 ■ = X 3 x + E ^ J s h(x 3 ,y a ,(i), 

s=0 

fc-1 

»=3 



(5.11) 



where j = 0,1, ... ,k. The operator $ is defined on the set 

Z(S) = {z=[(x j ,y j )% , \\z\\<8} , 
where the norm || • || is given as the maximum of modulus of components Xj, yj of the 



vector z. Notice that if zq = [(^,y?)]|_ is a fixed point of $, then the following 



diagram takes place 

(xg >V g) A... ^(xl,yl), 

i.e. the fixed point of $ gives a segment of an orbit of To. 

It is known [1] that, for small enough S = So and \x \ < So/2, \y k \ < So/2, the 
operator $ maps the set Z{5o) into itself and it is contracting. Thus, map (|5.11[) has 
a unique fixed point zq = [(x® (xq, yk), y® (xq, y k )]j=o that is limit of iterations under 
$ for any initial point from Z(6q). Thus, the coordinates x® and y° can be found by 
applying successive approximations. As an initial approximation, we take the solution 
of the linear problem: 

0(1) 0(1) j-k 

Xj = X J x , yf = 7 Vk 
It follows from (|5.11[) that the second approximation has a form 

xf 2 '' = \j Xo + E A^'- s - 1 A 2s 7 s - fc a;gyfeX 

((3 1 +O(\xf\x \ + \ 1 r k \y k \)) = 
Xx + AV fc E X-^'Y^Vkifa + O(\X\ s \x \ + \j\ s ~ k \y k \)) = 
A^o + (j - l^V'A- 1 ^ (/?! + O(|A| s |x | + \j\ s ~ k \y k \)) , 

yf 2) = i J - k y k + eV^-^V^zo^x 

(fa + O{\\\'\x \ + h\'- k \vk\)) = 
i J - k y k + i ] - 2k Ei-'^YxoyUfo + 0(\X\ s \x a \ + \i\ s - k \y k \)) = 

l 3 ~ k y k + (k- ih^-'xoyKh + O{\X\ s \x \ + h\ s ~ k \y k \)) 
Since 7 = A -1 , it follows from the precedent expression that 

\xf ) -Xixo\<L 1 jXi+\ 
\ y °W-X k -jy k \<L 2 (k-j)X 2k -j, 
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where L\ and L 2 are some positive constant independent of j and k. Substituting 
(|5.12p into (|5.11l) as the initial approximation, then the following ones will also satisfy 
estimates (|5.12[) . with the same constants Li and L 2 . Thus, formula (|3.5I) is valid for 
the coordinates x° and yg, fixed point of $. 

The estimates for the derivatives of the functions x® and j/g are deduced in the 
same way as done in [9] (see also modified versions of the proof in [38| [20" 1 121 ) ). 

6. Proofs of Lemmas [3], [4] and [5} 

6. 1 . Proof of Lemma [3| 

Since coordinates (£oi)2/oi) on a^ 1 are uniquely determined via cross-coordinates 
(^oi)2/n) in equations Q3.6p . we can express Tfc m as a map defined on points (xoi, j/n) 
and acting by the rule (xoi,yn) i— > (5oi)27ii)- As a result of this, we can express the 
map Tfe m in the following form 

xm - #2 = aA^oi + %n - yf) + '02(2/11 - yf) 2 + 
^ife(^oi,yn - y^jA*)) 

A™2/i2 (l + m\^hl l (y 12l x 02 ,fi)) =fi + c\ k x Q1 + d(y n - yf ) 2 + 
/iiA*a;oi(j/n - yf) + /o3(yn - yf) 3 + <f2k(x i, yn - 2/r>M)> 

X 2 n x 2 (l + m\ 2 n hl l (x 2,yi2, M)) = M + cA^yn + d(as i - yf) 2 + 
/11 Ail/ii (xoi - yf) + /o3(^oi - y^) 3 + <^2fe(yii - yf ^01, /^), 

yi2 - x$ = a ^iVn + b{xoi - Vi) + '02(^01 - yf) 2 + 
<^ifc(yn - yr,x i,M), 

where the coordinates xg 2 and j/12 are "intermediate" and 

<pi k {u,v,n) = O (\\ k u 2 + |AJ|H + M 3 ) , , . 

<P2k(u,v,») = 0(X\ k {u 2 + \u\v 2 ) + \X k \\u\v 2 ) +o(v 3 ). {K> - L > 

Now we perform the following shift in the coordinates 

£1 = xoi - xf + v\ m , 7}i = yn - yf + v\ m , 

6 = £02 - 3=2 + ^fem> »?1 = J/12 - 2/2~ + V lrm 

where vl — O(A^), i = 1,2, are some small coefficients which does not destroy the 
reversibility due to the condition (I3.5j) . Then, for suitable ft™, map Tfe m becomes 

6 = aAfCi + ^i + 'o2»7i + £ifc(£i,?7i>M), 
A 2 "r 72 (l + mA™^(r ?2 ,6,/i)) = 

A 2 ™6(l + mA 2 "^(6,r 7 2^)) = _ _ ""^ 

/x + cA^i + dCi/nA^iCi + /03^1 + <P2k(m> £l> 1"); 

?y 2 = aA^i + 6fi + ? 2^ 2 + Pifc(f?i,£i,M)> 

where, since relation (|3.5I) holds, we have that fi = [i— A 2 l (o! 2 + . . .) + cAi {a\ + . . .) and 
the new functions (p\k and (f>2k satisfy again conditions (|6.1I) . One must, however, to 
consider coefficients a, b, /g3 in (|6.2[) to be shifted by values of order O(fcAi) when 
comparing them with the initial coefficients in (|3.7I) . Substituting into the second and 
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third equations of (|6.2[) the expressions for £ 2 and V2 given by the first and the fourth 
precedent equations, we get an expression for T km of form 

XT (aXlm + + l 02 g) = 

p, + cAf£i + drj{ + JnA^i77i + /0377? + <P2k, 

XT (aA&i + 67?! + Z^f) = 

/2 + cAf?7i + d(( + /iiA^i^i + /03C1 + <P2k, 
which can be rewritten as 

aXTm + 6A' 2 "Ar fe 6 + io2AZ*Ar fc £? = 

AAr fc + cfi + dA- fe r;2 + / n ^ r/1 + f 03 X^ k r)l + X- k <p 2k , 

aXT^i + bXTX^ k m + i 02 XTK k Vi = U '" ' 

fiX^ k + <% + dXi k g + /um6 + /o 3 Ar fc Ci + Ar fe ^ 2 fe. 

Notice that the functions (p 2 k here may be changed in comparison with those in (|6.2p 
but still fulfill relations (|6.1[) . Finally, rescaling coordinates, 

6 = Xix, 771 = \\y, 

system (|6.3[) takes the form (|3.11l) where the coefficients c, d, . . . , Z02 are "original" 
ones (i.e., those appearing in formula (|3.7[1 ). 

6.2. Proof of Lemma^ 

The rescaled form (|3.11[) of the first-return map T km is, of course, implicit one and it 
corresponds to a formal representation Tkm '■ F(x, y) = G(x, y) which can be written 
in the explicit form (x, y) = Tk m (x, y) = F~ 1 G(x, y). Then we can find the Jacobian 
of T/cm using the relation 

D{T km )\ {x ^ D(F-')\ G{x y) D(G(x,y)), (6.4) 

where D(-) is the corresponding (differential) Jacobi matrix. At the fixed point 
(x = £*, y = 77*) of Tkm we can rewrite (16.41) as follows 

DG\ 



We find from (|3~TTj) that DF and DG are of the form 

/ 2dC + f llV *X k + 3s 03 X k (e) 2 c + / U £*A} 
^ bX^Xt? + ZXfloiC aXT 

( aXT bX^Xf + 2XTI02V* 

V c + /n^ A i' Zd^ + f^CXl+SXlsosiv*) 2 

plus terms of order 0(kX\ k ). Now we compute the Jacobian as 
J(T fcm )L, .,= detiDF-'DG) = d ^[ DG ] = 

-bcX^XT + 2adXTC - bfuXTC - 2d 02 A 2 "e + o(A^) (6 ' 5) 
-bcX^XT + 2adX k XTv* - bfn^fv* - 2cl a2 X 2 n V* + o(X k ) ' 
When the relation (|3.10[) is fulfilled we can rewrite (|6.5j) as 
T(T s, _ -bc + QX k r + o(A|) 

that gives relation f|3.20|) . 
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6.3. Proof of Lemma\^ 

Due to the reversibility, we can prove Lemma [5] directly for the truncated map 

H : x — M + cx — y 2 7 y — —-M + —y + —x 2 . We will use the following facts for 

c c c 

this map: it can be written in the explicit form (|2.2p and that for M > — j(c — l) 2 
it has a pair of symmetric fixed points P + = (p+,p+) and P~ — for which 

coordinates the formula (|3.18|) holds. Denote by p either p+ or p_ and let us assume 
that the corresponding fixed point P (i.e. P + or P~) is elliptic. Then, c and M 
have to take values from the open regions in the (c, M)-space of parameters given in 
Figure [5] 

The first step in our process is to shift the new origin of coordinates into the 
point (p, p) and to perform (Jordan) linear normal form, which leads our map to the 
following form 

. 2pcosip 2 

x = cos ip ■ x — sm ip ■ y — — — y + 

csm ip 

1 — 4p 2 — c cos ip i _ . , n - ,\ ^ 2\ 2 

^ a -9 9 ■ — ; (—csmip ■ x + (1 — ccosipw + 2py ) , 

4c p sm ip 



(6.6) 



2 P 2 

y = sm ip ■ x + cos ip ■ y y + 

1 - ° 2 

(— csixnp ■ x + (1 — ccosip)y + 2py 2 ) 



Acp 2 

where x, y, x, y stand again for the new variables. The linear part of (|6.6[) is a rotation 
of angle ip. At that point, we will lead our map into the so-called Birkhoff Normal 
Form up to order 3: z — e^z + d,2\z 2 z* + O4. To do it, we need to assume that 
A = e 1 ^ is not a Mh-root of unity for k = 3, 4 (the cases k = 1,2 correspond to 
parabolic fixed points and, thus, respond for boundaries of existence regions of elliptic 
fixed points). The coefficient B\ = —id2\e~ 1 ^' is called the first Birkhoff coefficient. 
By the Arnol'd-Moser Twist Theorem [39], the inequality B\ ^ (together with the 
absence of strong resonances) ensures that the elliptic point is generic or, in other 
words, KAM-stable. 

Introducing complex coordinates z — x + iy, z* = x — iy, map (|6.6|) takes the form 
z = e^z + A 20 z 2 + A llZ z* + A 02 (z*) 2 + A 21 z 2 z* + 4 (z, z*). 

Since we are assuming e 1 ^ not to be a 3rd or 4th root of unity (and also ip ^ 0,7r), 
our map can be lead into BNF up to order 3 and, afterwards, provides the following 
formula for the first Birkhoff coefficient 

(c + 1 - y (g + 1 + 2p)(c - 1 + 2p) 
325 4 psin 3 V'(2cosV' + 1) 

where 

Pi(c,p) = 64p 4 + 8(1 - 5)p 3 - 4(35 2 + Ac + 3)p 2 + 
2{c-l){c+l) 2 p-{c-l) 2 {c + l) 2 . 

Using this formul jptl we represent in Figure [9] curves B\{c,M) — 0, where the 
elliptic fixed point can be, a priori, not KAM-stable. In this Figure curves related to 
the strong resonances are presented. For resonances 1:1 and 1:2 the equations of the 

ft Notice that in the particular case c = —1, map H corresponds to H 2 , where H : x = y, y = 
M — x — y 2 is the Henon map. In this case we have that (i/>) = 2B^^ where ip = 2<p. It is not 
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Figure 9. In the (c, M)-plane, three grey and one hatching regions correspond to 
the existence of elliptic points: P + for the grey regions and P~ for the hatching 
one. Lines corresponding to the main resonances and vanishing the first Birkhoff 
coefficient for the elliptic point are shown and labelled. 



corresponding curves are given in Section I3T41 The equations of the 1:3 and 1:4 curves 
are as follows: 



c 2 + 1 v / c 2_ +T 
M = — - — ± (1 — c) for 1:4 resonance 



and 



iu g2 + 5 + 1 -u ^ + ~ c + 1 n ~\ t iq 

M = ± 1 — c for 1:3 resonance. 

4 2 K ' 



This completes the proof. 
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